Allosteric PI3Kα Inhibition Overcomes On-target Resistance to Orthosteric Inhibitors Mediated by Secondary PIK3CA Mutations.

PIK3CA mutations occur in ∼8% of cancers, including ∼40% of HR-positive breast cancers, where the PI3K-alpha (PI3Kα)-selective inhibitor alpelisib is FDA approved in combination with fulvestrant. Although prior studies have identified resistance mechanisms, such as PTEN loss, clinically acquired resistance to PI3Kα inhibitors remains poorly understood. Through serial liquid biopsies and rapid autopsies in 39 patients with advanced breast cancer developing acquired resistance to PI3Kα inhibitors, we observe that 50% of patients acquire genomic alterations within the PI3K pathway, including PTEN loss and activating AKT1 mutations. Notably, although secondary PIK3CA mutations were previously reported to increase sensitivity to PI3Kα inhibitors, we identified emergent secondary resistance mutations in PIK3CA that alter the inhibitor binding pocket. Some mutations had differential effects on PI3Kα-selective versus pan-PI3K inhibitors, but resistance induced by all mutations could be overcome by the novel allosteric pan-mutant-selective PI3Kα-inhibitor RLY-2608. Together, these findings provide insights to guide strategies to overcome resistance in PIK3CA-mutated cancers. SIGNIFICANCE: In one of the largest patient cohorts analyzed to date, this study defines the clinical landscape of acquired resistance to PI3Kα inhibitors. Genomic alterations within the PI3K pathway represent a major mode of resistance and identify a novel class of secondary PIK3CA resistance mutations that can be overcome by an allosteric PI3Kα inhibitor. See related commentary by Gong and Vanhaesebroeck, p. 204 . See related article by Varkaris et al., p. 240 . This article is featured in Selected Articles from This Issue, p. 201.

Latent human herpesvirus 6 is reactivated in CAR T cells.

Cell therapies have yielded durable clinical benefits for patients with cancer, but the risks associated with the development of therapies from manipulated human cells are understudied. For example, we lack a comprehensive understanding of the mechanisms of toxicities observed in patients receiving T cell therapies, including recent reports of encephalitis caused by reactivation of human herpesvirus 6 (HHV-6)1. Here, through petabase-scale viral genomics mining, we examine the landscape of human latent viral reactivation and demonstrate that HHV-6B can become reactivated in cultures of human CD4+ T cells. Using single-cell sequencing, we identify a rare population of HHV-6 'super-expressors' (about 1 in 300-10,000 cells) that possess high viral transcriptional activity, among research-grade allogeneic chimeric antigen receptor (CAR) T cells. By analysing single-cell sequencing data from patients receiving cell therapy products that are approved by the US Food and Drug Administration2 or are in clinical studies3-5, we identify the presence of HHV-6-super-expressor CAR T cells in patients in vivo. Together, the findings of our study demonstrate the utility of comprehensive genomics analyses in implicating cell therapy products as a potential source contributing to the lytic HHV-6 infection that has been reported in clinical trials1,6-8 and may influence the design and production of autologous and allogeneic cell therapies.

Germline mutations and developmental mosaicism underlying EGFR-mutant lung cancer.

While the development of multiple primary tumors in smokers with lung cancer can be attributed to carcinogen-induced field cancerization, the occurrence of multiple primary tumors in individuals with EGFR-mutant lung cancer who lack known environmental exposures remains unexplained. We identified ten patients with early-stage, resectable non-small cell lung cancer who presented with multiple anatomically distinct EGFR-mutant tumors. We analyzed the phylogenetic relationships among multiple tumors from each patient using whole exome sequencing (WES) and hypermutable poly-guanine (poly-G) repeat genotyping, as orthogonal methods for lineage tracing. In two patients, we identified germline EGFR variants, which confer moderately enhanced signaling when modeled in vitro. In four other patients, developmental mosaicism is supported by the poly-G lineage tracing and WES, indicating a common non-germline cell-of-origin. Thus, developmental mosaicism and germline variants define two distinct mechanisms of genetic predisposition to multiple EGFR-mutant primary tumors, with implications for understanding their etiology and clinical management.

Genetic events associated with venetoclax resistance in CLL identified by whole-exome sequencing of patient samples.

Although BCL2 mutations are reported as later occurring events leading to venetoclax resistance, many other mechanisms of progression have been reported though remain poorly understood. Here, we analyze longitudinal tumor samples from 11 patients with disease progression while receiving venetoclax to characterize the clonal evolution of resistance. All patients tested showed increased in vitro resistance to venetoclax at the posttreatment time point. We found the previously described acquired BCL2-G101V mutation in only 4 of 11 patients, with 2 patients showing a very low variant allele fraction (0.03%-4.68%). Whole-exome sequencing revealed acquired loss(8p) in 4 of 11 patients, of which 2 patients also had gain (1q21.2-21.3) in the same cells affecting the MCL1 gene. In vitro experiments showed that CLL cells from the 4 patients with loss(8p) were more resistant to venetoclax than cells from those without it, with the cells from 2 patients also carrying gain (1q21.2-21.3) showing increased sensitivity to MCL1 inhibition. Progression samples with gain (1q21.2-21.3) were more susceptible to the combination of MCL1 inhibitor and venetoclax. Differential gene expression analysis comparing bulk RNA sequencing data from pretreatment and progression time points of all patients showed upregulation of proliferation, B-cell receptor (BCR), and NF-κB gene sets including MAPK genes. Cells from progression time points demonstrated upregulation of surface immunoglobulin M and higher pERK levels compared with those from the preprogression time point, suggesting an upregulation of BCR signaling that activates the MAPK pathway. Overall, our data suggest several mechanisms of acquired resistance to venetoclax in CLL that could pave the way for rationally designed combination treatments for patients with venetoclax-resistant CLL.

Clinical Importance of the lncRNA NEAT1 in Cancer Patients Treated with Immune Checkpoint Inhibitors.

PURPOSE: mAbs targeting the PD-1/PD-L1 immune checkpoint are powerful tools to improve the survival of patients with cancer. Understanding the molecular basis of clinical response to these treatments is critical to identify patients who can benefit from this immunotherapy. In this study, we investigated long noncoding RNA (lncRNA) expression in patients with cancer treated with anti-PD-1/PD-L1 immunotherapy. EXPERIMENTAL DESIGN: lncRNA expression profile was analyzed in one cohort of patients with melanoma and two independent cohorts of patients with glioblastoma (GBM) undergoing anti-PD-1/PD-L1 immunotherapy. Single-cell RNA-sequencing analyses were performed to evaluate lncRNA expression in tumor cells and tumor-infiltrating immune cells. RESULTS: We identified the lncRNA NEAT1 as commonly upregulated between patients with melanoma with complete therapeutic response and patients with GBM with longer survival following anti-PD-1/PD-L1 treatment. Gene set enrichment analyses revealed that NEAT1 expression was strongly associated with the IFNγ pathways, along with downregulation of cell-cycle-related genes. Single-cell RNA-sequencing analyses revealed NEAT1 expression across multiple cell types within the GBM microenvironment, including tumor cells, macrophages, and T cells. High NEAT1 expression levels in tumor cells correlated with increased infiltrating macrophages and microglia. In these tumor-infiltrating myeloid cells, we found that NEAT1 expression was linked to enrichment in TNFα/NFκB signaling pathway genes. Silencing NEAT1 suppressed M1 macrophage polarization and reduced the expression of TNFα and other inflammatory cytokines. CONCLUSIONS: These findings suggest an association between NEAT1 expression and patient response to anti-PD-1/PD-L1 therapy in melanoma and GBM and have important implications for the role of lncRNAs in the tumor microenvironment.

Lesion Shedding Model: unraveling site-specific contributions to ctDNA.

Sampling circulating tumor DNA (ctDNA) using liquid biopsies offers clinically important benefits for monitoring cancer progression. A single ctDNA sample represents a mixture of shed tumor DNA from all known and unknown lesions within a patient. Although shedding levels have been suggested to hold the key to identifying targetable lesions and uncovering treatment resistance mechanisms, the amount of DNA shed by any one specific lesion is still not well characterized. We designed the Lesion Shedding Model (LSM) to order lesions from the strongest to the poorest shedding for a given patient. By characterizing the lesion-specific ctDNA shedding levels, we can better understand the mechanisms of shedding and more accurately interpret ctDNA assays to improve their clinical impact. We verified the accuracy of the LSM under controlled conditions using a simulation approach as well as testing the model on three cancer patients. The LSM obtained an accurate partial order of the lesions according to their assigned shedding levels in simulations and its accuracy in identifying the top shedding lesion was not significantly impacted by number of lesions. Applying LSM to three cancer patients, we found that indeed there were lesions that consistently shed more than others into the patients' blood. In two of the patients, the top shedding lesion was one of the only clinically progressing lesions at the time of biopsy suggesting a connection between high ctDNA shedding and clinical progression. The LSM provides a much needed framework with which to understand ctDNA shedding and to accelerate discovery of ctDNA biomarkers. The LSM source code has been available in the IBM BioMedSciAI Github (https://github.com/BiomedSciAI/Geno4SD).

Pirtobrutinib targets BTK C481S in ibrutinib-resistant CLL but second-site BTK mutations lead to resistance.

Covalent inhibitors of Bruton tyrosine kinase (BTK) have transformed the therapy of chronic lymphocytic leukemia (CLL), but continuous therapy has been complicated by the development of resistance. The most common resistance mechanism in patients whose disease progresses on covalent BTK inhibitors (BTKis) is a mutation in the BTK 481 cysteine residue to which the inhibitors bind covalently. Pirtobrutinib is a highly selective, noncovalent BTKi with substantial clinical activity in patients whose disease has progressed on covalent BTKi, regardless of BTK mutation status. Using in vitro ibrutinib-resistant models and cells from patients with CLL, we show that pirtobrutinib potently inhibits BTK-mediated functions including B-cell receptor (BCR) signaling, cell viability, and CCL3/CCL4 chemokine production in both BTK wild-type and C481S mutant CLL cells. We demonstrate that primary CLL cells from responding patients on the pirtobrutinib trial show reduced BCR signaling, cell survival, and CCL3/CCL4 chemokine secretion. At time of progression, these primary CLL cells show increasing resistance to pirtobrutinib in signaling inhibition, cell viability, and cytokine production. We employed longitudinal whole-exome sequencing on 2 patients whose disease progressed on pirtobrutinib and identified selection of alternative-site BTK mutations, providing clinical evidence that secondary BTK mutations lead to resistance to noncovalent BTKis.

Distinct cellular dynamics associated with response to CAR-T therapy for refractory B cell lymphoma.

Chimeric antigen receptor (CAR)-T cell therapy has revolutionized the treatment of hematologic malignancies. Approximately half of patients with refractory large B cell lymphomas achieve durable responses from CD19-targeting CAR-T treatment; however, failure mechanisms are identified in only a fraction of cases. To gain new insights into the basis of clinical response, we performed single-cell transcriptome sequencing of 105 pretreatment and post-treatment peripheral blood mononuclear cell samples, and infusion products collected from 32 individuals with large B cell lymphoma treated with either of two CD19 CAR-T products: axicabtagene ciloleucel (axi-cel) or tisagenlecleucel (tisa-cel). Expansion of proliferative memory-like CD8 clones was a hallmark of tisa-cel response, whereas axi-cel responders displayed more heterogeneous populations. Elevations in CAR-T regulatory cells among nonresponders to axi-cel were detected, and these populations were capable of suppressing conventional CAR-T cell expansion and driving late relapses in an in vivo model. Our analyses reveal the temporal dynamics of effective responses to CAR-T therapy, the distinct molecular phenotypes of CAR-T cells with differing designs, and the capacity for even small increases in CAR-T regulatory cells to drive relapse.

Molecular map of chronic lymphocytic leukemia and its impact on outcome.

Recent advances in cancer characterization have consistently revealed marked heterogeneity, impeding the completion of integrated molecular and clinical maps for each malignancy. Here, we focus on chronic lymphocytic leukemia (CLL), a B cell neoplasm with variable natural history that is conventionally categorized into two subtypes distinguished by extent of somatic mutations in the heavy-chain variable region of immunoglobulin genes (IGHV). To build the 'CLL map,' we integrated genomic, transcriptomic and epigenomic data from 1,148 patients. We identified 202 candidate genetic drivers of CLL (109 new) and refined the characterization of IGHV subtypes, which revealed distinct genomic landscapes and leukemogenic trajectories. Discovery of new gene expression subtypes further subcategorized this neoplasm and proved to be independent prognostic factors. Clinical outcomes were associated with a combination of genetic, epigenetic and gene expression features, further advancing our prognostic paradigm. Overall, this work reveals fresh insights into CLL oncogenesis and prognostication.

AXL and Error-Prone DNA Replication Confer Drug Resistance and Offer Strategies to Treat EGFR-Mutant Lung Cancer.

Anticancer therapies have been limited by the emergence of mutations and other adaptations. In bacteria, antibiotics activate the SOS response, which mobilizes error-prone factors that allow for continuous replication at the cost of mutagenesis. We investigated whether the treatment of lung cancer with EGFR inhibitors (EGFRi) similarly engages hypermutators. In cycling drug-tolerant persister (DTP) cells and in EGFRi-treated patients presenting residual disease, we observed upregulation of GAS6, whereas ablation of GAS6's receptor, AXL, eradicated resistance. Reciprocally, AXL overexpression enhanced DTP survival and accelerated the emergence of T790M, an EGFR mutation typical to resistant cells. Mechanistically, AXL induces low-fidelity DNA polymerases and activates their organizer, RAD18, by promoting neddylation. Metabolomics uncovered another hypermutator, AXL-driven activation of MYC, and increased purine synthesis that is unbalanced by pyrimidines. Aligning anti-AXL combination treatments with the transition from DTPs to resistant cells cured patient-derived xenografts. Hence, similar to bacteria, tumors tolerate therapy by engaging pharmacologically targetable endogenous mutators. SIGNIFICANCE: EGFR-mutant lung cancers treated with kinase inhibitors often evolve resistance due to secondary mutations. We report that in similarity to the bacterial SOS response stimulated by antibiotics, endogenous mutators are activated in drug-treated cells, and this heralds tolerance. Blocking the process prevented resistance in xenograft models, which offers new treatment strategies. This article is highlighted in the In This Issue feature, p. 2483.

Integrative molecular and clinical profiling of acral melanoma links focal amplification of 22q11.21 to metastasis.

Acral melanoma, the most common melanoma subtype among non-White individuals, is associated with poor prognosis. However, its key molecular drivers remain obscure. Here, we perform integrative genomic and clinical profiling of acral melanomas from 104 patients treated in North America (n = 37) or China (n = 67). We find that recurrent, late-arising focal amplifications of cytoband 22q11.21 are a leading determinant of inferior survival, strongly associated with metastasis, and linked to downregulation of immunomodulatory genes associated with response to immune checkpoint blockade. Unexpectedly, LZTR1 - a known tumor suppressor in other cancers - is a key candidate oncogene in this cytoband. Silencing of LZTR1 in melanoma cell lines causes apoptotic cell death independent of major hotspot mutations or melanoma subtypes. Conversely, overexpression of LZTR1 in normal human melanocytes initiates processes associated with metastasis, including anchorage-independent growth, formation of spheroids, and an increase in MAPK and SRC activities. Our results provide insights into the etiology of acral melanoma and implicate LZTR1 as a key tumor promoter and therapeutic target.

Paralog knockout profiling identifies DUSP4 and DUSP6 as a digenic dependence in MAPK pathway-driven cancers.

Although single-gene perturbation screens have revealed a number of new targets, vulnerabilities specific to frequently altered drivers have not been uncovered. An important question is whether the compensatory relationship between functionally redundant genes masks potential therapeutic targets in single-gene perturbation studies. To identify digenic dependencies, we developed a CRISPR paralog targeting library to investigate the viability effects of disrupting 3,284 genes, 5,065 paralog pairs and 815 paralog families. We identified that dual inactivation of DUSP4 and DUSP6 selectively impairs growth in NRAS and BRAF mutant cells through the hyperactivation of MAPK signaling. Furthermore, cells resistant to MAPK pathway therapeutics become cross-sensitized to DUSP4 and DUSP6 perturbations such that the mechanisms of resistance to the inhibitors reinforce this mechanism of vulnerability. Together, multigene perturbation technologies unveil previously unrecognized digenic vulnerabilities that may be leveraged as new therapeutic targets in cancer.

A common methodological phylogenomics framework for intra-patient heteroplasmies to infer SARS-CoV-2 sublineages and tumor clones.

BACKGROUND: All diseases containing genetic material undergo genetic evolution and give rise to heterogeneity including cancer and infection. Although these illnesses are biologically very different, the ability for phylogenetic retrodiction based on the genomic reads is common between them and thus tree-based principles and assumptions are shared. Just as the different frequencies of tumor genomic variants presupposes the existence of multiple tumor clones and provides a handle to computationally infer them, we postulate that the different variant frequencies in viral reads offers the means to infer multiple co-infecting sublineages. RESULTS: We present a common methodological framework to infer the phylogenomics from genomic data, be it reads of SARS-CoV-2 of multiple COVID-19 patients or bulk DNAseq of the tumor of a cancer patient. We describe the Concerti computational framework for inferring phylogenies in each of the two scenarios.To demonstrate the accuracy of the method, we reproduce some known results in both scenarios. We also make some additional discoveries. CONCLUSIONS: Concerti successfully extracts and integrates information from multi-point samples, enabling the discovery of clinically plausible phylogenetic trees that capture the heterogeneity known to exist both spatially and temporally. These models can have direct therapeutic implications by highlighting "birth" of clones that may harbor resistance mechanisms to treatment, "death" of subclones with drug targets, and acquisition of functionally pertinent mutations in clones that may have seemed clinically irrelevant. Specifically in this paper we uncover new potential parallel mutations in the evolution of the SARS-CoV-2 virus. In the context of cancer, we identify new clones harboring resistant mutations to therapy.

K-Mer Analyses Reveal Different Evolutionary Histories of Alpha, Beta, and Gamma Papillomaviruses.

Papillomaviruses (PVs) are a heterogeneous group of DNA viruses that can infect fish, birds, reptiles, and mammals. PVs infecting humans (HPVs) phylogenetically cluster into five genera (Alpha-, Beta-, Gamma-, Mu- and Nu-PV), with differences in tissue tropism and carcinogenicity. The evolutionary features associated with the divergence of Papillomaviridae are not well understood. Using a combination of k-mer distributions, genetic metrics, and phylogenetic algorithms, we sought to evaluate the characteristics and differences of Alpha-, Beta- and Gamma-PVs constituting the majority of HPV genomes. A total of 640 PVs including 442 HPV types, 27 non-human primate PV types, and 171 non-primate animal PV types were evaluated. Our analyses revealed the highest genetic diversity amongst Gamma-PVs compared to the Alpha and Beta PVs, suggesting reduced selective pressures on Gamma-PVs. Using a sequence alignment-free trimer (k = 3) phylogeny algorithm, we reconstructed a phylogeny that grouped most HPV types into a monophyletic clade that was further split into three branches similar to alignment-based classifications. Interestingly, a subset of low-risk Alpha HPVs (the species Alpha-2, 3, 4, and 14) split from other HPVs and were clustered with non-human primate PVs. Surprisingly, the trimer-constructed phylogeny grouped the Gamma-6 species types originally isolated from the cervicovaginal region with the main Alpha-HPV clade. These data indicate that characterization of papillomavirus heterogeneity via orthogonal approaches reveals novel insights into the biological understanding of HPV genomes.

Parallel Genomic Alterations of Antigen and Payload Targets Mediate Polyclonal Acquired Clinical Resistance to Sacituzumab Govitecan in Triple-Negative Breast Cancer.

Sacituzumab govitecan (SG), the first antibody-drug conjugate (ADC) approved for triple-negative breast cancer, incorporates the anti-TROP2 antibody hRS7 conjugated to a topoisomerase-1 (TOP1) inhibitor payload. We sought to identify mechanisms of SG resistance through RNA and whole-exome sequencing of pretreatment and postprogression specimens. One patient exhibiting de novo progression lacked TROP2 expression, in contrast to robust TROP2 expression and focal genomic amplification of TACSTD2/TROP2 observed in a patient with a deep, prolonged response to SG. Analysis of acquired genomic resistance in this case revealed one phylogenetic branch harboring a canonical TOP1 E418K resistance mutation and subsequent frameshift TOP1 mutation, whereas a distinct branch exhibited a novel TACSTD2/TROP2 T256R missense mutation. Reconstitution experiments demonstrated that TROP2T256R confers SG resistance via defective plasma membrane localization and reduced cell-surface binding by hRS7. These findings highlight parallel genomic alterations in both antibody and payload targets associated with resistance to SG. SIGNIFICANCE: These findings underscore TROP2 as a response determinant and reveal acquired SG resistance mechanisms involving the direct antibody and drug payload targets in distinct metastatic subclones of an individual patient. This study highlights the specificity of SG and illustrates how such mechanisms will inform therapeutic strategies to overcome ADC resistance.This article is highlighted in the In This Issue feature, p. 2355.

Longitudinal Single-Cell Dynamics of Chromatin Accessibility and Mitochondrial Mutations in Chronic Lymphocytic Leukemia Mirror Disease History.

While cancers evolve during disease progression and in response to therapy, temporal dynamics remain difficult to study in humans due to the lack of consistent barcodes marking individual clones in vivo. We employ mitochondrial single-cell assay for transposase-accessible chromatin with sequencing to profile 163,279 cells from 9 patients with chronic lymphocytic leukemia (CLL) collected across disease course and utilize mitochondrial DNA (mtDNA) mutations as natural genetic markers of cancer clones. We observe stable propagation of mtDNA mutations over years in the absence of strong selective pressure, indicating clonal persistence, but dramatic changes following tight bottlenecks, including disease transformation and relapse posttherapy, paralleled by acquisition of copy-number variants and changes in chromatin accessibility and gene expression. Furthermore, we link CLL subclones to distinct chromatin states, providing insight into nongenetic sources of relapse. mtDNA mutations thus mirror disease history and provide naturally occurring genetic barcodes to enable patient-specific study of cancer subclonal dynamics. SIGNIFICANCE: Single-cell multi-omic profiling of CLL reveals the utility of somatic mtDNA mutations as in vivo barcodes, which mark subclones that can evolve over time along with changes in accessible chromatin and gene expression profiles to capture dynamics of disease evolution. See related commentary by Hilton and Scott, p. 2965. This article is highlighted in the In This Issue feature, p. 2945.

Explainable AI reveals changes in skin microbiome composition linked to phenotypic differences.

Alterations in the human microbiome have been observed in a variety of conditions such as asthma, gingivitis, dermatitis and cancer, and much remains to be learned about the links between the microbiome and human health. The fusion of artificial intelligence with rich microbiome datasets can offer an improved understanding of the microbiome's role in human health. To gain actionable insights it is essential to consider both the predictive power and the transparency of the models by providing explanations for the predictions. We combine the collection of leg skin microbiome samples from two healthy cohorts of women with the application of an explainable artificial intelligence (EAI) approach that provides accurate predictions of phenotypes with explanations. The explanations are expressed in terms of variations in the relative abundance of key microbes that drive the predictions. We predict skin hydration, subject's age, pre/post-menopausal status and smoking status from the leg skin microbiome. The changes in microbial composition linked to skin hydration can accelerate the development of personalized treatments for healthy skin, while those associated with age may offer insights into the skin aging process. The leg microbiome signatures associated with smoking and menopausal status are consistent with previous findings from oral/respiratory tract microbiomes and vaginal/gut microbiomes respectively. This suggests that easily accessible microbiome samples could be used to investigate health-related phenotypes, offering potential for non-invasive diagnosis and condition monitoring. Our EAI approach sets the stage for new work focused on understanding the complex relationships between microbial communities and phenotypes. Our approach can be applied to predict any condition from microbiome samples and has the potential to accelerate the development of microbiome-based personalized therapeutics and non-invasive diagnostics.

The Lipogenic Regulator SREBP2 Induces Transferrin in Circulating Melanoma Cells and Suppresses Ferroptosis.

Circulating tumor cells (CTC) are shed by cancer into the bloodstream, where a viable subset overcomes oxidative stress to initiate metastasis. We show that single CTCs from patients with melanoma coordinately upregulate lipogenesis and iron homeostasis pathways. These are correlated with both intrinsic and acquired resistance to BRAF inhibitors across clonal cultures of BRAF-mutant CTCs. The lipogenesis regulator SREBP2 directly induces transcription of the iron carrier Transferrin (TF), reducing intracellular iron pools, reactive oxygen species, and lipid peroxidation, thereby conferring resistance to inducers of ferroptosis. Knockdown of endogenous TF impairs tumor formation by melanoma CTCs, and their tumorigenic defects are partially rescued by the lipophilic antioxidants ferrostatin-1 and vitamin E. In a prospective melanoma cohort, presence of CTCs with high lipogenic and iron metabolic RNA signatures is correlated with adverse clinical outcome, irrespective of treatment regimen. Thus, SREBP2-driven iron homeostatic pathways contribute to cancer progression, drug resistance, and metastasis. SIGNIFICANCE: Through single-cell analysis of primary and cultured melanoma CTCs, we have uncovered intrinsic cancer cell heterogeneity within lipogenic and iron homeostatic pathways that modulates resistance to BRAF inhibitors and to ferroptosis inducers. Activation of these pathways within CTCs is correlated with adverse clinical outcome, pointing to therapeutic opportunities.This article is highlighted in the In This Issue feature, p. 521.